A binary trait model reveals the fitness effects of HIV-1 escape from T cell responses.

Natural selection often acts on multiple traits simultaneously. For example, the virus HIV-1 faces pressure to evade host immunity while also preserving replicative fitness. While past work has studied selection during HIV-1 evolution, it is challenging to quantitatively separate different contributions to fitness. This task is made more difficult because a single mutation can affect both immune escape and replication. Here, we develop an evolutionary model that disentangles the effects of escaping CD8+ T cell-mediated immunity, which we model as a binary trait, from other contributions to fitness. After validation in simulations, we applied this model to study within-host HIV-1 evolution in a clinical data set. We observed strong selection for immune escape, sometimes greatly exceeding past estimates, especially early in infection. Conservative estimates suggest that roughly half of HIV-1 fitness gains during the first months to years of infection can be attributed to T cell escape. Our approach is not limited to HIV-1 or viruses, and could be adapted to study the evolution of quantitative traits in other contexts.

Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data.

In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.

popDMS infers mutation effects from deep mutational scanning data.

Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions.

An inference model gives insights into innate immune adaptation and repertoire diversity.

The innate immune system is the body's first line of defense against infection. Natural killer (NK) cells, a vital part of the innate immune system, help to control infection and eliminate cancer. Studies have identified a vast array of receptors that NK cells use to discriminate between healthy and unhealthy cells. However, at present, it is difficult to explain how NK cells will respond to novel stimuli in different environments. In addition, the expression of different receptors on individual NK cells is highly stochastic, but the reason for these variegated expression patterns is unclear. Here, we studied the recognition of unhealthy target cells as an inference problem, where NK cells must distinguish between healthy targets with normal variability in ligand expression and ones that are clear "outliers." Our mathematical model fits well with experimental data, including NK cells' adaptation to changing environments and responses to different target cells. Furthermore, we find that stochastic, "sparse" receptor expression profiles are best able to detect a variety of possible threats, in agreement with experimental studies of the NK cell repertoire. While our study was specifically motivated by NK cells, our model is general and could also apply more broadly to explain principles of target recognition for other immune cell types.

Small RNA sequencing of field Culex mosquitoes identifies patterns of viral infection and the mosquito immune response.

Mosquito-borne disease remains a significant burden on global health. In the United States, the major threat posed by mosquitoes is transmission of arboviruses, including West Nile virus by mosquitoes of the Culex genus. Virus metagenomic analysis of mosquito small RNA using deep sequencing and advanced bioinformatic tools enables the rapid detection of viruses and other infecting organisms, both pathogenic and non-pathogenic to humans, without any precedent knowledge. In this study, we sequenced small RNA samples from over 60 pools of Culex mosquitoes from two major areas of Southern California from 2017 to 2019 to elucidate the virome and immune responses of Culex. Our results demonstrated that small RNAs not only allowed the detection of viruses but also revealed distinct patterns of viral infection based on location, Culex species, and time. We also identified miRNAs that are most likely involved in Culex immune responses to viruses and Wolbachia bacteria, and show the utility of using small RNA to detect antiviral immune pathways including piRNAs against some pathogens. Collectively, these findings show that deep sequencing of small RNA can be used for virus discovery and surveillance. One could also conceive that such work could be accomplished in various locations across the world and over time to better understand patterns of mosquito infection and immune response to many vector-borne diseases in field samples.

Bézier interpolation improves the inference of dynamical models from data.

Many dynamical systems, from quantum many-body systems to evolving populations to financial markets, are described by stochastic processes. Parameters characterizing such processes can often be inferred using information integrated over stochastic paths. However, estimating time-integrated quantities from real data with limited time resolution is challenging. Here, we propose a framework for accurately estimating time-integrated quantities using Bézier interpolation. We applied our approach to two dynamical inference problems: Determining fitness parameters for evolving populations and inferring forces driving Ornstein-Uhlenbeck processes. We found that Bézier interpolation reduces the estimation bias for both dynamical inference problems. This improvement was especially noticeable for data sets with limited time resolution. Our method could be broadly applied to improve accuracy for other dynamical inference problems using finitely sampled data.

Estimating linkage disequilibrium and selection from allele frequency trajectories.

Genetic sequences collected over time provide an exciting opportunity to study natural selection. In such studies, it is important to account for linkage disequilibrium to accurately measure selection and to distinguish between selection and other effects that can cause changes in allele frequencies, such as genetic hitchhiking or clonal interference. However, most high-throughput sequencing methods cannot directly measure linkage due to short-read lengths. Here we develop a simple method to estimate linkage disequilibrium from time-series allele frequencies. This reconstructed linkage information can then be combined with other inference methods to infer the fitness effects of individual mutations. Simulations show that our approach reliably outperforms inference that ignores linkage disequilibrium and, with sufficient sampling, performs similarly to inference using the true linkage information. We also introduce two regularization methods derived from random matrix theory that help to preserve its performance under limited sampling effects. Overall, our method enables the use of linkage-aware inference methods even for data sets where only allele frequency time series are available.

Inferring Epistasis from Genetic Time-series Data.

Epistasis refers to fitness or functional effects of mutations that depend on the sequence background in which these mutations arise. Epistasis is prevalent in nature, including populations of viruses, bacteria, and cancers, and can contribute to the evolution of drug resistance and immune escape. However, it is difficult to directly estimate epistatic effects from sampled observations of a population. At present, there are very few methods that can disentangle the effects of selection (including epistasis), mutation, recombination, genetic drift, and genetic linkage in evolving populations. Here we develop a method to infer epistasis, along with the fitness effects of individual mutations, from observed evolutionary histories. Simulations show that we can accurately infer pairwise epistatic interactions provided that there is sufficient genetic diversity in the data. Our method also allows us to identify which fitness parameters can be reliably inferred from a particular data set and which ones are unidentifiable. Our approach therefore allows for the inference of more complex models of selection from time-series genetic data, while also quantifying uncertainty in the inferred parameters.

Pathogenesis, Symptomatology, and Transmission of SARS-CoV-2 through Analysis of Viral Genomics and Structure.

The novel coronavirus SARS-CoV-2, which emerged in late 2019, has since spread around the world and infected hundreds of millions of people with coronavirus disease 2019 (COVID-19). While this viral species was unknown prior to January 2020, its similarity to other coronaviruses that infect humans has allowed for rapid insight into the mechanisms that it uses to infect human hosts, as well as the ways in which the human immune system can respond. Here, we contextualize SARS-CoV-2 among other coronaviruses and identify what is known and what can be inferred about its behavior once inside a human host. Because the genomic content of coronaviruses, which specifies the virus's structure, is highly conserved, early genomic analysis provided a significant head start in predicting viral pathogenesis and in understanding potential differences among variants. The pathogenesis of the virus offers insights into symptomatology, transmission, and individual susceptibility. Additionally, prior research into interactions between the human immune system and coronaviruses has identified how these viruses can evade the immune system's protective mechanisms. We also explore systems-level research into the regulatory and proteomic effects of SARS-CoV-2 infection and the immune response. Understanding the structure and behavior of the virus serves to contextualize the many facets of the COVID-19 pandemic and can influence efforts to control the virus and treat the disease. IMPORTANCE COVID-19 involves a number of organ systems and can present with a wide range of symptoms. From how the virus infects cells to how it spreads between people, the available research suggests that these patterns are very similar to those seen in the closely related viruses SARS-CoV-1 and possibly Middle East respiratory syndrome-related CoV (MERS-CoV). Understanding the pathogenesis of the SARS-CoV-2 virus also contextualizes how the different biological systems affected by COVID-19 connect. Exploring the structure, phylogeny, and pathogenesis of the virus therefore helps to guide interpretation of the broader impacts of the virus on the human body and on human populations. For this reason, an in-depth exploration of viral mechanisms is critical to a robust understanding of SARS-CoV-2 and, potentially, future emergent human CoVs (HCoVs).

Pathogenesis, Symptomatology, and Transmission of SARS-CoV-2 through Analysis of Viral Genomics and Structure.

The novel coronavirus SARS-CoV-2, which emerged in late 2019, has since spread around the world and infected hundreds of millions of people with coronavirus disease 2019 (COVID-19). While this viral species was unknown prior to January 2020, its similarity to other coronaviruses that infect humans has allowed for rapid insight into the mechanisms that it uses to infect human hosts, as well as the ways in which the human immune system can respond. Here, we contextualize SARS-CoV-2 among other coronaviruses and identify what is known and what can be inferred about its behavior once inside a human host. Because the genomic content of coronaviruses, which specifies the virus's structure, is highly conserved, early genomic analysis provided a significant head start in predicting viral pathogenesis and in understanding potential differences among variants. The pathogenesis of the virus offers insights into symptomatology, transmission, and individual susceptibility. Additionally, prior research into interactions between the human immune system and coronaviruses has identified how these viruses can evade the immune system's protective mechanisms. We also explore systems-level research into the regulatory and proteomic effects of SARS-CoV-2 infection and the immune response. Understanding the structure and behavior of the virus serves to contextualize the many facets of the COVID-19 pandemic and can influence efforts to control the virus and treat the disease.

Adenovirus-vectored vaccine containing multidimensionally conserved parts of the HIV proteome is immunogenic in rhesus macaques.

An effective vaccine that can protect against HIV infection does not exist. A major reason why a vaccine is not available is the high mutability of the virus, which enables it to evolve mutations that can evade human immune responses. This challenge is exacerbated by the ability of the virus to evolve compensatory mutations that can partially restore the fitness cost of immune-evading mutations. Based on the fitness landscapes of HIV proteins that account for the effects of coupled mutations, we designed a single long peptide immunogen comprising parts of the HIV proteome wherein mutations are likely to be deleterious regardless of the sequence of the rest of the viral protein. This immunogen was then stably expressed in adenovirus vectors that are currently in clinical development. Macaques immunized with these vaccine constructs exhibited T-cell responses that were comparable in magnitude to animals immunized with adenovirus vectors with whole HIV protein inserts. Moreover, the T-cell responses in immunized macaques strongly targeted regions contained in our immunogen. These results suggest that further studies aimed toward using our vaccine construct for HIV prophylaxis and cure are warranted.

MPL resolves genetic linkage in fitness inference from complex evolutionary histories.

Genetic linkage causes the fate of new mutations in a population to be contingent on the genetic background on which they appear. This makes it challenging to identify how individual mutations affect fitness. To overcome this challenge, we developed marginal path likelihood (MPL), a method to infer selection from evolutionary histories that resolves genetic linkage. Validation on real and simulated data sets shows that MPL is fast and accurate, outperforming existing inference approaches. We found that resolving linkage is crucial for accurately quantifying selection in complex evolving populations, which we demonstrate through a quantitative analysis of intrahost HIV-1 evolution using multiple patient data sets. Linkage effects generated by variants that sweep rapidly through the population are particularly strong, extending far across the genome. Taken together, our results argue for the importance of resolving linkage in studies of natural selection.

Predominance of positive epistasis among drug resistance-associated mutations in HIV-1 protease.

Drug-resistant mutations often have deleterious impacts on replication fitness, posing a fitness cost that can only be overcome by compensatory mutations. However, the role of fitness cost in the evolution of drug resistance has often been overlooked in clinical studies or in vitro selection experiments, as these observations only capture the outcome of drug selection. In this study, we systematically profile the fitness landscape of resistance-associated sites in HIV-1 protease using deep mutational scanning. We construct a mutant library covering combinations of mutations at 11 sites in HIV-1 protease, all of which are associated with resistance to protease inhibitors in clinic. Using deep sequencing, we quantify the fitness of thousands of HIV-1 protease mutants after multiple cycles of replication in human T cells. Although the majority of resistance-associated mutations have deleterious effects on viral replication, we find that epistasis among resistance-associated mutations is predominantly positive. Furthermore, our fitness data are consistent with genetic interactions inferred directly from HIV sequence data of patients. Fitness valleys formed by strong positive epistasis reduce the likelihood of reversal of drug resistance mutations. Overall, our results support the view that strong compensatory effects are involved in the emergence of clinically observed resistance mutations and provide insights to understanding fitness barriers in the evolution and reversion of drug resistance.

Cross-serotypically conserved epitope recommendations for a universal T cell-based dengue vaccine.

Dengue virus (DENV)-associated disease is a growing threat to public health across the globe. Co-circulating as four different serotypes, DENV poses a unique challenge for vaccine design as immunity to one serotype predisposes a person to severe and potentially lethal disease upon infection from other serotypes. Recent experimental studies suggest that an effective vaccine against DENV should elicit a strong T cell response against all serotypes, which could be achieved by directing T cell responses toward cross-serotypically conserved epitopes while avoiding serotype-specific ones. Here, we used experimentally-determined DENV T cell epitopes and patient-derived DENV sequences to assess the cross-serotypic variability of the epitopes. We reveal a distinct near-binary pattern of epitope conservation across serotypes for a large number of DENV epitopes. Based on the conservation profile, we identify a set of 55 epitopes that are highly conserved in at least 3 serotypes. Most of the highly conserved epitopes lie in functionally important regions of DENV non-structural proteins. By considering the global distribution of human leukocyte antigen (HLA) alleles associated with these DENV epitopes, we identify a potentially robust subset of HLA class I and class II restricted epitopes that can serve as targets for a universal T cell-based vaccine against DENV while covering ~99% of the global population.

Inference of compressed Potts graphical models.

We consider the problem of inferring a graphical Potts model on a population of variables. This inverse Potts problem generally involves the inference of a large number of parameters, often larger than the number of available data, and, hence, requires the introduction of regularization. We study here a double regularization scheme, in which the number of Potts states (colors) available to each variable is reduced and interaction networks are made sparse. To achieve the color compression, only Potts states with large empirical frequency (exceeding some threshold) are explicitly modeled on each site, while the others are grouped into a single state. We benchmark the performances of this mixed regularization approach, with two inference algorithms, adaptive cluster expansion (ACE) and pseudolikelihood maximization (PLM), on synthetic data obtained by sampling disordered Potts models on Erdos-Rényi random graphs. We show in particular that color compression does not affect the quality of reconstruction of the parameters corresponding to high-frequency symbols, while drastically reducing the number of the other parameters and thus the computational time. Our procedure is also applied to multisequence alignments of protein families, with similar results.

Deconvolving mutational patterns of poliovirus outbreaks reveals its intrinsic fitness landscape.

Vaccination has essentially eradicated poliovirus. Yet, its mutation rate is higher than that of viruses like HIV, for which no effective vaccine exists. To investigate this, we infer a fitness model for the poliovirus viral protein 1 (vp1), which successfully predicts in vitro fitness measurements. This is achieved by first developing a probabilistic model for the prevalence of vp1 sequences that enables us to isolate and remove data that are subject to strong vaccine-derived biases. The intrinsic fitness constraints derived for vp1, a capsid protein subject to antibody responses, are compared with those of analogous HIV proteins. We find that vp1 evolution is subject to tighter constraints, limiting its ability to evade vaccine-induced immune responses. Our analysis also indicates that circulating poliovirus strains in unimmunized populations serve as a reservoir that can seed outbreaks in spatio-temporally localized sub-optimally immunized populations.

MPF-BML: a standalone GUI-based package for maximum entropy model inference.

SUMMARY: Learning underlying correlation patterns in data is a central problem across scientific fields. Maximum entropy models present an important class of statistical approaches for addressing this problem. However, accurately and efficiently inferring model parameters are a major challenge, particularly for modern high-dimensional applications such as in biology, for which the number of parameters is enormous. Previously, we developed a statistical method, minimum probability flow-Boltzmann Machine Learning (MPF-BML), for performing fast and accurate inference of maximum entropy model parameters, which was applied to genetic sequence data to estimate the fitness landscape for the surface proteins of human immunodeficiency virus and hepatitis C virus. To facilitate seamless use of MPF-BML and encourage more widespread application to data in diverse fields, we present a standalone cross-platform package of MPF-BML which features an easy-to-use graphical user interface. The package only requires the input data (protein sequence data or data of multiple configurations of a complex system with large number of variables) and returns the maximum entropy model parameters. AVAILABILITY AND IMPLEMENTATION: The MPF-BML software is publicly available under the MIT License at https://github.com/ahmedaq/MPF-BML-GUI. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Modelling and in vitro testing of the HIV-1 Nef fitness landscape.

An effective vaccine is urgently required to curb the HIV-1 epidemic. We have previously described an approach to model the fitness landscape of several HIV-1 proteins, and have validated the results against experimental and clinical data. The fitness landscape may be used to identify mutation patterns harmful to virus viability, and consequently inform the design of immunogens that can target such regions for immunological control. Here we apply such an analysis and complementary experiments to HIV-1 Nef, a multifunctional protein which plays a key role in HIV-1 pathogenesis. We measured Nef-driven replication capacities as well as Nef-mediated CD4 and HLA-I down-modulation capacities of thirty-two different Nef mutants, and tested model predictions against these results. Furthermore, we evaluated the models using 448 patient-derived Nef sequences for which several Nef activities were previously measured. Model predictions correlated significantly with Nef-driven replication and CD4 down-modulation capacities, but not HLA-I down-modulation capacities, of the various Nef mutants. Similarly, in our analysis of patient-derived Nef sequences, CD4 down-modulation capacity correlated the most significantly with model predictions, suggesting that of the tested Nef functions, this is the most important in vivo. Overall, our results highlight how the fitness landscape inferred from patient-derived sequences captures, at least in part, the in vivo functional effects of mutations to Nef. However, the correlation between predictions of the fitness landscape and measured parameters of Nef function is not as accurate as the correlation observed in past studies for other proteins. This may be because of the additional complexity associated with inferring the cost of mutations on the diverse functions of Nef.

Characterization of Intact Proviruses in Blood and Lymph Node from HIV-Infected Individuals Undergoing Analytical Treatment Interruption.

The role of lymphoid tissue as a potential source of HIV-1 rebound following interruption of antiretroviral therapy (ART) is uncertain. To address this issue, we compared the latent viruses obtained from CD4+ T cells in peripheral blood and lymph nodes to viruses emerging during treatment interruption. Latent viruses were characterized by sequencing near-full-length (NFL) proviral DNA and env from viral outgrowth assays (VOAs). Five HIV-1-infected individuals on ART were studied, four of whom participated in a clinical trial of a TLR9 agonist that included an analytical treatment interruption. We found that 98% of intact or replication-competent clonal sequences overlapped between blood and lymph node. In contrast, there was no overlap between 205 latent reservoir and 125 rebound sequences in the four individuals who underwent treatment interruption. However, rebound viruses could be accounted for by recombination. The data suggest that CD4+ T cells carrying latent viruses circulate between blood and lymphoid tissues in individuals on ART and support the idea that recombination may play a role in the emergence of rebound viremia.IMPORTANCE HIV-1 persists as a latent infection in CD4+ T cells that can be found in lymphoid tissues in infected individuals during ART. However, the importance of this tissue reservoir and its contribution to viral rebound upon ART interruption are not clear. In this study, we sought to compare latent HIV-1 from blood and lymph node CD4+ T cells from five HIV-1-infected individuals. Further, we analyzed the contribution of lymph node viruses to viral rebound. We observed that the frequencies of intact proviruses were the same in blood and lymph node. Moreover, expanded clones of T cells bearing identical proviruses were found in blood and lymph node. These latent reservoir sequences did not appear to be the direct origin of rebound virus. Instead, latent proviruses were found to contribute to the rebound compartment by recombination.